Foraging habitat and performances of King penguins 
Aptenodytes patagonicus Miller, 1778 at Kerguelen islands 
in relation to climatic variability 

by 

Charles-Andre BOST* * (1), Anne GOARANT (2), Annette SCHEFFER (1), 

Philippe KOUBBI (2), Guy DUHAMEF (3) & Jean-Benoit CHARRASSIN (3) 


ABSTRACT. - Understanding how climate change affects the foraging ecology of key marine predators is an important 
issue in the study of Southern Ocean food webs. Since 1998, we have conducted a long term research program on the 
foraging ecology of a top-diving predator, the King penguin Aptenodytes patagonicus, simultaneously at Kerguelen and 
Crozet Islands with the support of the IPEV (Institut Polaire Franjais). King penguin is one of the most important predatory 
birds at Kerguelen. In summer, they forage preferentially along large-scale physical features and use the three dimensions 
of hydrographic features to feed on myctophid fishes (the penguins’ main prey and a major prey species in the Southern 
Ocean). Primarily, this study aims to evaluate how the change in foraging parameters and success of penguins reflect the 
impact of oceanographic conditions on key food webs in the polar frontal zone. Each summer, breeding (incubation and 
brooding) penguins from the Ratmanoff colony (Courbet Peninsula, Kerguelen; 100 000 pairs) are instrumented with Argos 
transmitters or GPS with time-temperature-depth recorders. Here we provide a first modelling approach to the King pen¬ 
guin’s foraging habitat during summer from static variables (bathymetry), dynamic variables (SSHt, SST, Chlorophyll) and 
their related trends (gradients). In addition, trawl data on the distribution of myctophid fishes have been integrated. The 
most important factors explaining penguin foraging location are the SSHt, the bathymetry gradient and the SST. Results of 
interannual change in penguins at-sea trajectories, diving behaviour, foraging success and effort are discussed according to 
the interannual changes in the hydrographic structure in the Kerguelen region. 


RESUME. - Habitat d’alimentation du manchot royal ( Aptenodytes patagonicus) a Kerguelen. 

Comprendre comment les changements climatiques vont affecter les predateurs superieurs marins de l’ocean Austral 
et les chames trophiques dont ils dependent est d’un interet majeur. Depuis 1998, nous menons un programme de recher¬ 
che a long terme sur les strategies alimentaires d’un predateur marin plongeur, le manchot royal Aptenodytes patagonicus 
aux lies Kerguelen et Crozet, avec le support de 1’IPEV (Institut Polaire Frangais). Le manchot royal est un des predateurs 
marins les plus importants de Kerguelen, en termes de biomasse consommante. En ete, les manchots royaux se nourrissent 
preferentiellement de poissons de la famille des myctophides, un groupe cle des chames trophiques de l’ocean Austral. Ils 
recherchent leurs proies preferentiellement au niveau de caracteristiques physiques a grande echelle et utilisent les 3 dimen¬ 
sions de la structure hydrologique. Un des objectifs de notre etude est d’evaluer comment les changements dans l’ecologie 
alimentaire des manchots vont refleter l’impact des conditions oceanographiques a l’echelle locale et a grande echelle, 
sur les chames alimentaires dont ils dependent. A cette fin, des manchots partant en mer pour se nourrir sont equipes de 
balises Argos, GPS ou enregistreurs de plongee / temperature de l’eau, durant chaque ete austral. A Kerguelen, la colonie 
etudiee est situee a Ratmanofff, sur la fagade est de la peninsule Courbet (100 000 couples). Nous presentons ici la premiere 
approche de modelisation de l’habitat alimentaire en mer des manchots, a partir de l’analyse des variables statiques (bathy- 
metrie) et dynamiques (temperature et hauteur de l’eau et leurs gradients, chlorophylle et gradients associes). Les donnees 
des peches scientifiques existantes sur la distribution des myctophides ont ete integrees. Les facteurs explicatifs les plus 
importants vis-a -vis de la distribution en mer des manchots sont les anomalies des hauteurs d’eau et les gradients associes, 
la temperature de l’eau et le gradient de bathymetrie. Les resultats des changements d’annee en annee des trajectoires en 
mer des manchots, du comportement de plongee et de leur effort de peche sont discutes en perspectives des changements 
hydrologiques affectant l’ecosysteme de Kerguelen. 
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The Southern Ocean is one of the marine areas most sus¬ 
ceptible to climate change (Trathan et al., 2007). The detec¬ 
tion and quantification of the impact of the climatic vari¬ 


ability on the biodiversity of this ocean and its food webs 
have been given high international priority because it still 
accommodates the planet’s largest unexploited animal bio- 
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mass (Atkinson et al ., 2004). However, the links between the 
physical features of the Southern Ocean, biological produc¬ 
tivity and the distribution and abundance of zooplanktonic 
and nektonic prey remain poorly understood. 

Considering these difficulties, it has been proposed that 
the foraging behaviour of tagged predators could be used 
as bio-indicators of the availability of underlying resource 
distribution (Wilson et al., 1994; Bost et al., 1997; Wilson 
et al., 2002; Austin et al., 2006). Many marine birds from 
these areas are wide-ranging predators, highly mobile, and 
dependent on secondary and tertiary productivities. At sea, 
they strive to forage efficiently to maximize their chances of 
reproductive success (Lescroel et al., 2010). 

Most research efforts concerning pelagic resources in the 
Southern Oceans have been devoted to the change in the sta¬ 
tus of krill stocks (Atkinson et al., 2004). However, the mes- 
opelagic fishes (Family Myctophidae) are also an important 
biological resource, in which biomass is estimated at 200- 
400 x 10 6 tons (e.g., Lubimova et al., 1987; Pakhomov etal., 
1994). Additionally the distribution of these small schooling 
fish is closely related to the thermal structure of the water 
mass (Torres and Somero, 1988). Myctophids are difficult to 
sample using traditional techniques because of their patchy 
distribution and mobility (Duhamel 1998). Diving top pred¬ 
ators such as King penguins Aptenodytes patagonicus Miller 
1778 are good candidates to investigate the inter-annual and 
long-term change in myctophid distribution. King penguins 
are one of the most important avian consumers of the South¬ 
ern Ocean (Guinet et al., 1996). It is also one of the most 
specialized seabirds in terms of diet, relying almost exclu¬ 
sively on myctophid fishes during the summer (Cherel and 
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Figure 1. - Foraging dives distribution of King penguins satellite- 
tracked from Ratmanoff Colony, East of Kerguelen Islands (n = 4 
years, 1998-2002,27 271 dives, n = 27 birds). 


Ridoux, 1992). They are deep divers and pelagic foragers 
able to routinely dive deeper than 250 m and to forage up to 
400 km from their colonies (Bost et al., 2002). 

Since 1998, we have conducted a long term research 
program on the King penguin foraging ecology at Kergue¬ 
len Islands. The aim of the project is to evaluate how the 
changes in foraging parameters and success of an avian top 
predator can reflect and predict change in the availability of 
myctophids in one of the most productive ecosystem of the 
Polar Frontal Zone, the Kerguelen archipelago. 

Here we provide a first modelling approach of the King 
penguins foraging habitat during summer at Kerguelen. The 
habitat model was developed to explain and predict spatial 
distribution of foraging effort within the penguins’ available 
geographic range. This was carried out from the analysis 
of penguins foraging activity {via a bio-logging approach) 
over four years (1999-2002) concurrently with both datasets 
describing physical and biological oceanography. 


MATERIALS AND METHODS 

The project relied on the long-term monitoring of pen¬ 
guins movements at sea and foraging effort during the breed¬ 
ing season depending on oceanographic conditions. Each 
summer, breeding (incubation and brooding) penguins from 
the Ratmanoff colony (Courbet Peninsula, 100 000 pairs) are 
instrumented with Argos transmitters or Fast-loc GPS (Sir- 
track: Havelock North, NZ), and Time-Temperature-Depth 
recorders (MK7 to MK9, Wildlife Computers: Redmond, 
WA, USA). 

We used the spatial distribution of dives to determine the 
at-sea distribution of the foraging habitat. Only dives deeper 
than 50 m were used as most of the feeding activity occurs 
beyond that depth (Charrassin et al., 2002a). This corre¬ 
sponded to a total of 27 271 foraging dives recorded on 27 
birds from 1998 to 2002. A grid of 0.2° cell size was designed 
over the Kerguelen shelf and in each cell the number of dive 
was recorded. Only cells where at least three different indi¬ 
viduals had dived were kept (n = 107, 54% of the dives; Fig. 
1). The number of foraging dives was then log-transformed. 
Finally the number of dives per grid cell was modelled 
according to three types of environmental components: the 
spatial, physical and biological component. Dynamic ocea¬ 
nographic variables (surface height, sea surface temperature, 
chlorophyll concentration) were extracted for each grid cell 
from both sources of oceanographic information (model 
and satellite), whereas static variables such as bathymetry 
was obtained from the National Oceanic and Atmospheric 
Administration’s (NOAA) ETOPO dataset, and their gradi¬ 
ent (in terms of slope between each cell) computed. 

In addition, myctophid presence-absence (distribution 
trawls data on the Kerguelen shelf from the Ichtyoker data- 
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base, 1998-2000) was modelled according to the four pre¬ 
vious environmental variables and their gradient. A model 
was performed using Random Forests, a machine-learning 
method combining three models trained on different boot¬ 
strap replicate sample of the data (Elith and Graham, 2009). 
We used data from night trawls at depths ranging from the 
surface to 100 m (n=73 trawls), performed from January to 
March (Duhamel 1998; Guinet et al., 2001). The perform¬ 
ance and evaluation of the model was evaluated using a 
Leave-one-out cross-validation. 

The variable testing procedure for assessing penguin for¬ 
aging habitat was based on the test of the main variables and 
their possible combinations, which provided seven hypoth- 
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Figure 2. - Output of the foraging habitat model showing the most 
important variables (Hypothesis “Space + environment”) .VI: large 
scale distribution of dives (Kerguelen North and South); V4: mes- 
oscale distribution of dives (range 100 km, including the 2 areas of 
highest diving activity); V6: mesoscale distribution of dives (range 
76 km, limit of the studied area); V14: mesoscale distribution of 
dives (range 52 km, restricted diving zones, Kerguelen North); 
V26: mesoscale distribution of dives (range 47 km); V34: fine scale 
distribution of dives (range < 10 km) [not shown]. 


eses (Tab. I). A bootstrap procedure (Potts and Elith, 2006) 
was performed to evaluate the performance of each model¬ 
ling hypothesis (10 000 simulations). 

RESULTS AND DISCUSSION 

At Kerguelen, King penguins are distributed in an area 
dominated by complex interactions between the proximity 
of the polar front, bathymetry and current advection. They 
forage at a mean distance of 300 km off the colony and tar¬ 
get mainly waters along the shelf break that are associated 
with the Polar Front. We used the most important physi¬ 
cal and biological variables and 
parameters likely to determine the 
penguin’s foraging distribution. 

The three main influencing 
variables follow. First, the spa¬ 
tial distribution of dives, which 
are strongly autocorrelated at dif¬ 
ferent scales. To determine the 
corresponding spatial scales, a 
principal coordinate analysis of 
neighbour matrices (PCNM) was 
performed (Dray et al., 2006) 
(Tab. I). Six spatial scales were 
retained in the analysis (Fig. 2). 
Second, the environmental vari¬ 
ables, which includes the physi¬ 
cal (static: bathymetry; dynamic: 
SST, sea surface height anomaly) 
and biological variables (chloro¬ 
phyll). Lastly, the prey fields (myctophids) obtained from 
the model. 

Overall, the hypothesis combining spatial and environ¬ 
mental variables had the best average performance (45% of 
explained variance. Tab. I). The output of the foraging habi¬ 
tat model indicated the contribution of the most significant 
variables (Fig. 2), primarily, the variable combining the 
spatial distribution of dives and environmental variables; 
and secondly, two significant dynamic variables, sea surface 
height and sea surface temperature. Static variables, like 
depth gradient and depth, also play a significant role. 

Intrinsically, modelling of the King penguin’s foraging 
habitat based on a four-year study of diving activity indicates 
the importance of physical variables such as the temperature 
of the water mass and sea surface height. Other studies have 
shown how oceanic foragers such as King penguins concen¬ 
trate their foraging effort in colder waters of the Polar Fron¬ 
tal Zone where myctophids tend to aggregate (Cotte et al., 
2007; Bost et al., 1997). In the Crozet sector, the foraging 
activity occurs within mesoscale frontal zones and strong 
currents, both associated with eddies at the Polar Front (Cotte 


Table I. - Tests of the hypothesis predicting the distribution of the king penguin foraging habitat 
at Kerguelen Islands (RMSE: Root Means Square Error; AVE: Average Error). The letters m 
and b correspond to the coefficients of a linear regression between the observed and predicted 
values: observed = m x predicted + b. The hypothesis with the best performance is indicated in 
bold. 


Hypothesis 

Pear¬ 

son 

Spear¬ 

man 

m 

h 

RMSE 

AVE 

Variance 

explained 

“Space” 

0.652 

0.631 

1.232 

-0.467 

0.259 

0.000 

27.9% 

“Physical Environment” 

0.667 

0.487 

1.502 

0.025 

0.236 

0.004 

28.1% 

“Prey” 

0.481 

0.142 

1.017 

-1.326 

0.278 

-0.001 

17.3% 

“Space + Environment” 

0.787 

0.703 

1.217 

-0.448 

0.203 

0.004 

45.1% 

“Space + Prey” 

0.741 

0.647 

1.311 

-0.635 

0.222 

0.001 

39.3% 

“Environment + Prey” 

0.665 

0.487 

1.515 

-0.044 

0.238 

0.006 

29.1% 

“Space + environment + 
prey” 

0.784 

0.690 

1.251 

-0.518 

0.203 

0.004 

43.7% 
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et al., 2007). However, this study also confirms the role of 
the extensive Kerguelen plateau as a major physical vari¬ 
able driving the foraging behaviour of such pelagic forager 
(Charrassin et al., 2002b). Climate change scenario predicts 
a major southern shift of the polar front during this century 
(Solomon et al., 2007). In the future, warmer surface waters 
and increased frequency of mesoscale warm anomalies may 
alter the myctophid distribution and the foraging success of 
the penguins and other top-diving predators. Next steps will 
be to model the foraging habitat/prey distribution in distinct 
climatic scenarios and to evaluate the foraging responses of 
the penguins in different modelled foraging habitats by using 
Behaviour Based Models (Grimm and Railsback. 2005). 
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